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ABSTRACT 

The integrated luminosity of the TP-AGB phase is a major uncertainty in stel- 
lar population synthesis models. We use the white dwarf initial final mass relation 
and stellar interiors models to demonstrate that a significant fraction of the core mass 
growth for intermediate (1.5 < M & < 6) mass stars takes place during the TP-AGB 
phase. We find evidence that the peak fractional core mass contribution for TP-AGB 
stars is ~ 20% and occurs for stars between 2 M Q and 3.5 M & . Using a simple fuel 
consumption argument we couple this core mass increase to a lower limit on the TP- 
AGB phase energy output. Roughly half of the energy released in models of TP-AGB 
stars can be directly accounted for by this core growth; while the remainder is predom- 
inantly the stellar yield of He. A robust measurement of the emitted light in this phase 
will therefore set strong constraints on helium enrichment from TP-AGB stars, and we 
estimate the yields predicted by current models as a function of initial mass. Implica- 
tions for stellar population studies and prospects for improvements are discussed. 



1. Introduction 

Stars can experience a core-collapse supernova only if they are born with a mass much higher 
than the Sun, even though a chemically evolved core of order only one solar mass is required to 
ignite the advanced burning stages. The culprit is mass loss severe enough to strip off the stellar 
envelope before a sufficiently massive processed core develops. Although there is some mass loss 
even in shell hydrogen burning giants, the vast majority in intermediate mass stars occurs in the 
presence of thermal pulsations involving interactions between hydrogen and helium burning shells; 
we refer to this as the thermally pulsing AGB (TP-AGB) phase. Our understanding of even the 
basic properties of the TP-AGB phase, such as the lifetime or light emitted, is limited because their 
severe mass loss is difficult to constrain observationally or predict theoretically. 

The uncertainties in TP-AGB evolution have profound consequences for stellar population 
studies. Stellar interiors models can accurately predict evolutionary properties up to the onset of 
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the TP-AGB phase, and such models have been extensively used to reconstruct star formation his- 
tories in both resolved and unresolved populations. Models of the TP-AGB phase, however, are 
extremely dependent on the assumptions related to mass loss and their predictive power is thus 
limited. There has been intriguing work suggesting t hat the fraction of l ight emitted by TP-AGB 
stars in LMC star clusters is very high, of order 40 % (IPersson et al.ll 19831) . Such a substantial flux, 
largely redistributed to the fa r-IR by dust can b e especially important for interpreting intermediate 
redshift galaxy properties (IConroy et al.ll2009|) . There has been a strong emphasis on updatin g 
population synthesis t o include the TP-AGB phase (e.g. lBruzual & Charlotll2003l : lMarastonll2005l) . 
Maraston et al.l (120061) demonstrated that modeling of the TP-AGB phas e has become a de fining 
characteristic of different stellar population synthesis (SPS) codes while lConroy et al.1 (I2009|) iden- 
tified the substantial uncertainties in TP-AGB properties as a major component of the error budget 
for galaxy evolution models. 

In this paper we employ a fuel consumption argument to set a firm lower bound on the fraction 
of light emitted during the TP-AGB phase. We use stellar interiors models to set the core mass 
at the onset of the TP-AGB phase, and demonstrate that it is surprisingly insensitive to the choice 
of input physics. The nuclear processed core then grows until the envelope is expelled, at which 
point the final white dwarf mass is set. The white dwarf initial-final mass relationship (IFMR) can 
in turn be inferred from open clusters. We use the difference between the final and starting masses 
in the TP-AGB phase as a bound on the fuel consumed, and thus the emitted light. Although it 
is common for investigators to compare their models to initial-final mass relationships, and fuel 
consumption arguments have been used to check population synthesis models, the quantitative 
bounds from the IFMR are typically not used in population synthesis calculations. We argue that 
the additional empirical informa tion encoded there permits a narrower set of possibilities than 
considered by IConroy et al.l (|2009|) . Processed fuel (especially helium) can be ejected in winds, so 
the white dwarf data formally sets only a lower bound on the emitted light. This raises interesting 
links between the light emitted in this phase and chemical evolution studies, which we discuss in 
our conclusions. The plan of our paper is straightforward. Our sample and methods are discussed 
in Section 2, our results are presented in Section 3, and we discuss their broader implications in 
Section 4. 



2. Sample and Methods 

Stellar evolution theory makes a robust prediction for the core mass at the onset of the TP- 
AGB phase as a function of composition and initial mass. Open clusters provide a laboratory 
where we can constrain the initial mass given a sample of known stellar remnants. Through fuel 
consumption arguments, the difference between the final core mass (Mi), set by the mass of the 
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white dwarf, and the core mass at the onset of the TP-AGB phase (M c j T p) is a lower limit on the 
total light emitted from stars during their TP-AGB phase as a function of initial mass. We therefore 
begin by discussing the white dwarf IFMR, which provides the largest and most accurate sample 
of initial-final mass pairs. We then address the theoretical interior models and their uncertainties, 
and close with how we relate fuel consumption to energy output. 

Our sample consists of 48 white dwarfs (WDs) in 9 open cluster s. The primary reference for 



the majority of our (M r M f ) data is the compilation of iSalaris et all (|2009); hereafter S09. S09 
provide an extensive investigation of the systematic errors involved in determining the IFMR and 
calculate initial and final masses self-consistently across their entire sample, i.e., they use the same 
set of physics to compute the WD masses, cluster ages, and initial masses. We limit our analysis 
to clusters appearing in the S09 study. We include initial and final mass pairs from nine of the ten 
open clusters in S09. The Pleiades is omitted because the white dwarf in question is very young, 
the progenitor was very massive, and the theoretical errors in M c j T p are large. We have updated the 
S09 WD measurements or cluster parameters (such as age or metallicity) for several clusters in our 
sample; we discuss the adopted WD and cluster data in Section [2TTI and Section [231 respectively. 



2.1. White Dwarf Samples and Final Masses 



Final white dwarf masses are determined from spectra of their Balmer lines, which are partic- 
ularly sensitive to changes in surface gravities (g) and effective temperature (T e ff). Theoretical WD 
cooling curves provide a mass-radius relationship which is used to infer the final white dwarf mass 
and the cooling age (r coo i) of the remnant, with only a mild dependence on the WD composition. 



S09 use the cooling curves of ISalaris et alj (I2000|) to determine M t and r cool and their uncer- 
tainties from these data. Investigating the systematic uncertainties inherent in the determination of 
the IFMR was a major focus of S09. As such, their reported error in Mf incorporates not only the 
uncertainty from obs ervational errors but also the resulting range of Mf calculated using alternative 
cooling tracks from lAlthaus et al. ( 2003 ) and limiting case choices of: neutrino energy loss rates, 
conductive opacities, WD core composition, and hydrogen envelope thickness. The average frac- 
tional error in final mass and cooling age for the S09 sample is 8.4% and 34.9%, respectively. Note 
that these errors are extremely conservative and incorporate many more potential error sources than 
in typical IFMR studies. 

We adopt the WD masses and cooling ages from S09 for six of the nine clusters in our sample. 
For these clusters, S09 has compiled the most recent and high resolution data. S09 adopt the log g 
and r eff mea surements first compiled i n lFerrario et a l. (2005) fo r the Hyades (IClaver et al.ll200l|) . 
NGC 2516 dKoester & Reimerslll996l) . and M37 dKalirai et alJl2005h . WD measurements from 
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NGC 6819 and NGC 7789 ( Kalirai et all 120081) : and NGC1039 dRubin et all 120081) complete the 
list of clusters whose WD data is up to date in S09. 

There are more recent WD observations in NGC 3532, M35, and Praesepe. In NGC3532, 
Dobbie et al.l (|20091) present high resolution sp ectroscopic and pho t ometr ic observations of six 
WDs, including the three originally observed in iKoester & Reimersl (| 19931) and reported in S09. 
With precise distance modulus measurements, they determine that two of these six WDs are not 
associated with the cluster. The final masses and cooling times for the remaining four WDs are 
interpolated from the theoretical cooling curves of Fontaine et al.1 (|2001l) . The average fractional 
error in final mass and cooling age is 6.0% and 22.2%, respectively. 



Williams et al.l (|2009a|) expand upon their previous study ([Williams et al.1 120041 reported in 
S09) of WDs in M35. They present high resolution spectroscopic observations and updated log g 
and r eff measurements of 12 DA WDs. Three WDs in their sample were not fit satisfactorily 
by model atmospheres (po ssibly due to magnetic fields) and w e remove them from our sam- 
ple. IWilliams et al.l (l2009al) interpolate the WD cooling curves of Fontaine et al.l (|200ll) to obtain 
masses and cooling times for the WDs in their sample. The mean fractional error in final mass for 
the nine WDs we analyze is 9.3% and is 51% in cooling time. 

Casewell et al. ( 2009|) present high-resolution spectroscopic observations of 9 WDs in Prae- 



sepe using the Very Large Telescope. They determine that the object WD0836+201 has been mis- 
labeled by previous studies and has strong magnetic fields. Using radial velocity measurements, 
they argue against the inclusion of another candidate, WD0837+218, as a cluster member. We use 
their remaining seven candidates as our Praesepe WDs. Taking advantage of their high resolution, 
Casewell et al.l (|2009l) measured line core velocity shifts of the Ha and Hf3 lines and refit model 



atmo spheres to the spectra accounting for these shifts, resulting in lower x 2 fits. iCasewell et al. 
(200 9j) interpo l ate th eir measured r eff and log g amongst a grid of cooling curve models from 



Fontaine et al.l (1200 1|) to determine the WD masses and cooling times. The average fractional 
uncertainty is 5.6% in WD mass and 1 1.6% in cooling age. 

We note our adopted WD masses from these sources are calculated using a different set of 
theoretical cooling curves from the rest of our sample. However, the uncertainty in WD mass is 
dominated by ob servational errors rather than systematic ones stemming from the choice of WD 
cooling curves (|Salaris et al. 2009 ). On the other hand, WD cooling ages vary more when inter- 
polated using different models. To account for this syste matic uncertainty, we assume a fractional 
error of 50% for the WD cooling ages obtained using the [Fontaine et al] (|200l|) models. For nomi- 
nal choices of WD composition, envelope thickness, energy loss rates, and opacities, the fractional 
uncertainty in WD cooling ages for the remainder of our sample is still < 20% in most cases. 
Improvements in the observational measurements of log g and T e s and the removal of cluster non- 
members from the sample more than offset any systematic affect stemming from the inclusion of 
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a second group of WD cooling curve models in our analysis. 



2.2. Cluster Parameters 

Progenitor lifetimes are a function of WD cooling ages and cluster ages. The accuracy of our 
initial-final mass data set depends critically on constraining the cluster distance and age. Although 
a few clusters have measured parallaxes, the distances to most are inferred by main sequence fitting. 
Photometric methods can be used to infer composition, reddening, and distance even for clusters 
with limited membership and spectroscopic data (see for example Pinsonneault et al. 1997, 2004; 
An et al. 2007). Well-studied open clusters can have spectroscopic metallicity measurements, ex- 
tinction inferred from polarization studies, and both radial velocity and proper motion membership 
data. The uncertainties in these basic cluster parameters are largely determined by the available 
information in each specific cluster, which we discuss below, and set the uncertainty in the inferred 
cluster distance. The measured distance to the cluster yields the turnoff luminosity, which can be 
combined with isochrones to determine the age and the mass-main sequence lifetime relationship 
(see S09 for a discussion). The errors in these two aspects are distinct in nature. 



Core size deser ves special comment. Co nvective core overshooting (IChiosi & Maederl ll986) 



or rotational mixing (M aeder & Meynetll2000l) are difficult to model theoretically, and both mech- 



anisms have the practical effect of extending the main sequence lifetime by providing extra fuel. 
The main sequence is broader than that predicted by models without oyershqot or mixing, whic h 



is evidence that this phenomenon is real to some degree (e.g. lAndersenlll991l : lTorres et al.1 12010). 
As a result, the main sequence lifetime-mass relationship has substantial theoretical uncertainties. 
However, the total fuel burned is more reliable than the cluster ages because the helium core mass 
at the end of the main sequence is less than the minimum core mass required for helium ignition. 
Uncertainties related to overshooting will thus be important for intermediate mass stars, where 
the star leaves the main sequence with a helium core greater than the minimum required for core 
He burning, but the impact on stars at the lower end of the range will be reduced. For the pur- 
poses of this paper we adopt a conservative limiting case approach, independently inferring fuel 
consumption bounds from models with and without overshoot. 



2.3. Initial Masses 

The progenitor lifetime is the difference between the cluster age (r c i us ) and the cooling age 
of the WD (r coo i). Stellar evolutionary tracks spanning a range of stellar mass and metallicity are 
then interpolated to yield the best-fitting progenitor mass (Mi). WD spectra modeled by theoreti- 



-6- 



cal cooling curves constrain r coo i (Section [27Tb . Set by the distance to the cluster and metallicity, 
the cluster turn-off luminosity is a direct indicator of r clus . The cluster color-magnitude diagram 
(CMD), [Fe/H], and E(B—V) are the three main data inputs in main sequence fitting algorithms 
used to constrain the distance modulus to the cluster. S09 first assume, using measurements from 
the literature, a metallicity and reddening for each cluster (see their Table 2 and 3 for their sources 
and values, respectively). To obtain the cluster distance, they empirically fit the main sequence of 
V-(B - V) CMDs using a sample of of field dwarfs with known metallicities and Hipparcos paral- 
laxes (IPercival et al.ll2003|) . Once the distance is known, the turn-off luminosity i s measured. S09 



interp olate the isochrones (both with and without convective overshooting) of IPietrinferni et al 



(I2004|) in both turn-off luminosity and metallicity to constrain the best fit cluster age. 



We adopt the cluster parameters, namely cluster age and composition, presented in S09 for 
three of the nine clusters in our analysis. The composition measurements included in S09 of the 
Hyades, NGC 6819, and NGC 1039 are either taken from the most recent, high resolution stud- 
ies of these clusters or are in agreement with current measurements. Addi tionally, the computed 



distances to these clust ers align with previou s distance investigations (e.g. IPerryman et al.lll998 



Kalirai et al.ll200ll : and lJones & Prosserll 19961 for the Hyades, NGC 6819, and NGC 1039 respec- 
tively). 

We searched the literature for revised metallicity measurements of and/or distance determina- 
tions to all the clusters in our sample. For six of the clusters, we modify the cluster composition, 
distance, or both from the values found in S09. In two of these six clusters, S09 either had diffi- 
culty with their main sequence (MS) fitting technique due to poor data (in the case of NGC 3532) 
or chose a reddening value significantly different than that calculated in the recent literature (M37). 
In other clusters in our sample, several high resolution spectroscopic studies over the last decade 
have produced very precise metallicity measurements. We use these more recent values in our final 
calculations. 



We put all of our cluster ages on the same relative scale. IPercival et al.1 (120031) present the 
main sequence fitting technique applied in S09, and derive a metallicity dependency of A(B — 
V) = 0A54A[Fe/H] for their procedure. Their algorithm produces distances relative to an as- 
sumed distance to the Hyades, principally by matching the shape of the Hyades MS to the cluster 
CM Ds. Assuming a line ar slope with a magnitude of ~ 5 for the relevant portion of the Hyades 



MS dPercivaletalJl2003h . A(B-V) = AM v /5; therefore, AM,, = 0.77 * A[Fe/H]. iTwarog et al. 



i2009h performed an investigation of MS fitting techniques using nearby field stars with Hippar- 
cos parallaxes and p recise metallic i ties an d fi nd AM V = 0.98 * A\F e/H]. The metallicity depen- 
dency determ i ned by IPercival et al] (12003b and lTwarog et al] (|2009l) are similar; implementing the 
ITwarog et al. ( 2009 ) relation would not significantly impact our results. We modify cluster ages 
corresponding to changes in distance according to the errors presented in S09. We assume that, for 
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a given cluster, the fractional error in distance is equal to the fractional error in age. If A(m-M) 
is the S09 cluster distance subtracted from the new distance then A(m-M)/a^ m - M ) * cr tclm = At c i us , 
where cr (m _M) and a Tdus are from S09 and the new cluster age is t dus = t dus S09 - At c i us . We adopt ages 
relative to those of S09 for consistency. Below we summarize our revisions to these six clusters. 



Praesepe lAn et al.l (120081) present a new distance measurement to this cluster using empiri- 
cally calibrated isochrones. Additionally, they obtain high signal to noise spectroscopy of several 
Praesepe stars, reporting [Fe/HJ=JXll ±JX03. This is between the [Fe/H]= 0.04 ±0.06 from 
the literature compilat ion of Grattonl d20o"o) used by S09 and the subsequent higher [Fe/H]= 
0.27 ±0.10 found bv IPace et all J2008h . lAnetaD J2008h find (m-M) = 6.33 ±0.04 assum- 
ing [Fe/H]= 0.14 ±0.02 (a weighted mean of their result and litera t ure v alues excluding non- 
members) and E(B-V) = 0.006 ±0.002. We adopt the latter ! An et all (120071) value for our metal- 
licity. Comparing the metallicity used in this distance determination with that of S09, we find 
A[Fe/H] = 0. 14 - 0. 04 = 0.10. Accounti ng for this change in metallicity according to the metal- 
licit); '_de2£^^ncejn2M^^M^MluP03), the de-reddened S09 distance is (m-M)o = 6.32 ±0.04. 
The An et al.l ((2008) distance is 0.01 mag larger than expected, implying a relatively younger 
cluster age (r clus ) of 637 ± 50 Myr using isochrones incorporating convective overshoot (OS) and 
440 ± 40 Myr assuming no overshoot (nOS). The errors in the cluster age remain unchanged. 



NGC 2516 lAnet all (120081) simultaneously best-fit their photometry of NGC 2516 with [Fe/H]= 
-0.04±0.05, E(B-V) = 0.1 17 ±0.002, and (m-Af) = 8.03 ±0.04. This represents a substantial 
improvement in crpe/H] over S09 (cr [Fe / H] = 0.11). De-reddening and accounting for A[Fe/H] = 
0.12, the S09 distance is (m-M) = 8.04. Hence, A(m - M) = -0.01 mag yielding lifetimes of 
t s = 137 ± 29 Myr and r NO s = 91 ± 26 Myr. The error in cluster age has been reduced from their 
S09 values by the factor a (m - M )AnO&/^(m-M),so9 = 0.57. 

NGC 3532 S09 could not find sufficiently precise photometry of this cluster in the magnitude 
range necessary to match their WD templates; subsequently, they had to scale the derived distance 



and age from their results for Praesepe. Assuming the same reddening as S09, iKharchenko et al 
(|2005|) find (m-M)y = 8.61 ±0.2 using optical photometry and Hipparcos proper motions. Their 
distance measurement is 0.21 mag larger than S09, giving r os = 316 ± 80 Myr and t NO s = 216 ± 
80 Myr. Note that the errors in cluster age are reduced by (J( m -M),K05 /°"(m-M),so9 = 0.80. We could 
not find a a recent, high resolution s pectroscopic metallicity determination of this cluster. We use 
the S09 value: [Fe/H]= 0.02 ±0.06 (brattonlboool) . 



M37 S09 fit the main sequence of this cluster using two different sets of cluster parameters. 
The first, with super-solar metallicity and E(B-V) = 0.30 produced a cluster age that was im- 
probably young. The progenitor masses for this age were so large relative to the remnant masses 
that the theoretically predicted core mass at the onset of the TP-AGB phase was typically greater 
than the remnant mass. S09 tried a substantially lower metallicity, [Fe/H]= -0.20, and reddening, 
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E(B— V) = 0.23 and obtained a more reasonable cluster age in agreemen t with previou s age d eter- 
minations. We scale our results with this latter cluster characteristic set. iHartman et al.l (|2008l) per- 
formed high resolution spectroscopic and photometric observations of M37 using the MMT. They 
determine [Fe/H]=0.045 ±0.044, E(B - V) = 0.227 ±0.038, and (m-M) v = 11.57 ±0.13. Com- 
paring this distance with that of S09, A(m-M) = 0. 17. However, adjusting the S09 distance to the 
same metallicity (A[Fe/H] = 0.24), we find {m-M) v = 11.58. Thus, A(m-M) = 1 1.57 - 1 1.58 = 
-0.01 mag, and the cluster age for M37 is r os = 554 ± 54 Myr and t nos = 354 ± 43 Myr. The errors 
in cluster age have increased by the factor cr (m _ M ),A«08 / o~(m-M),so9 = 1 -083 over their S09 counterparts. 

M35 We found one high re solution, high signal to noise m easurement of M35's composi- 
tion ([Fe/H]= -0.21 ± 0. 10) by iBarrado v Navascues et all (boOll) . In M35, A[Fe/H] = -0.21 - 
(-0. 19) = -0.02 and thus A(m-M) = -0.015. The revised cluster ages are slightly older than those 
in S09: t os = 124 ± 30 Myr and t nos = 88 ± 25 Myr. Cluster age errors remain unchanged from 
the S09 values. 



NGC 7789 iTautvaisiene et all (120051) determine [Fe/H]= -0.04 ± 0.05 in NGC 7789 using 
high resolution spectra of evolved cluster members while the measurement ([Fe/H]= -0.13 ±0.08) 
quoted in SQ9 is an a verage of photometric and low dispersion spectroscopic abundance measure- 
ments (|Gratt on 2000). A[Fe/H] =-0.04 -(-0.1 3) = 0.09; the accompanying change in distance is 
A(m-M) = 0.07. The cluster age is 58 Myr younger than stated in S09: t os = 1442 ± 100 Myr and 
tnos = 1042 ± 100 Myr. Cluster age errors are repeated from S09. 

We note that the largest departure in cluster age from S09 is 21%. This maximum deviation 
occurs in the case of NGC 3532, which has more limited data in S0 9 than any other cluster in 
the sample. Our adopted age for NGC 3532 is consistent with that of iKoester & Reimersl (|1993L 
T s = 302 ±154 Myr). Typically, however, our cluster ages represent modifications to those of S09 
on the order of a few percent. To reduce our error budget, we used more recent cluster composi- 
tion determinations for many of these clusters. We have demonstrated that our choices of cluster 
metallicity do not seriously impact our derived cluster ages. Our results would not have changed 
signifigantly had we simply adopted all the reported cluster ages in S09. 

For each WD, the lifetime of its progenitor is r prog = r clus -r mr ,i. Errors are propag ated from 
r cool and r clus to r prog . Using the stellar evolutionary models of IPietrinferni et al.1 (I2004|) . we con- 
struct a data cube of Mj as a function of r prog and [Fe/H]. For the six clusters with modified cluster 
ages or WD cooling ages, we interpolate the cube linearly in [Fe/H] and quadratically in r piog 
to obtain the progenitor mass. For each object, we interpolate this grid 2 x 10 6 times, each time 
independently sampling a Gaussian distribution in M± and [Fe/H] (the width of each Gaussian is 
<7 T and u [Fe/H] 5 respectively). One sigma error bars encapsulate the central 68% of the values 
interpolated fro m the grid. The uncert ainty in M; incorporates the error in r prog and [Fe/H]. S09 
interpolate the IPietrinferni et al.1 (|2004l) models in the same fashion to obtain progenitor masses for 
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WDs in the other three clusters. In their error budget, they also include the systematic offset in M[ 
from interpolating a different stellar evolutionary model set (namely those of the PADOVA group 
(IGirardi et alJboOOl) ). While our revised M\ measurements in Praesepe, NGC 2516, NGC 3532, 
M37, M35, and NGC 7789 do not include the systematic error associated with using the PADOVA 
models, this potential contribution to the total a Ml is small when compared to the cluster parameter 
uncertainties (S09). 



2.4. Cluster Averages 

The final and initial masses of all 48 WDs in our sample are shown in Figure [T] In many of the 
clusters, there is considerable scatter in these values, especially when M ; > 3M . When r c i us -r cool 
is small, errors in either quantity contribute more strongly to T prog , resulting in a spread of initial 
masses for a single cluster. We red uce this scatter by representing each cluster by its weighted 



mean M\ and Mf (first proposed in iKalirai et al.l 120081) . To determine the influence of outlying 



measurements, we also compute the median M ; and M f for each cluster and adopt the standard 
error approximation a med = n/2 * a mecm . We find no statistically significant change in our results 
when adopting the median Mi, Mf for each cluster instead of the mean. For the remainder of this 
paper we present our results using the weighted mean M; and Mf of each cluster. 



2.5. Fraction of Core Growth in the TP-AGB Phase 

The core mass at the first thermal pulse of the AGB stage is fixed by theory via the star's 
initial mass and composition. Similar to the calculation of M ; from r prog , we con struct a data cube 



of the core mass at the first thermal pulse (M c i TP ) given M; and [Fe/H] using the [Pietrinferni et al 



(|2004|) models. M c i TP is interpolated linearly in [Fe/H] and quadratically in M x . Central values 
and errors are calculated similarly to that of M\ in Section [231 The error in M c iT p is the one sigma 
range of interpolated M c 1TP given the uncertainties in cluster [Fe/H] and M x . As the final remnant 
mass is equal to the core mass at the tip of the AGB, the core grows by AM c =M f -M c 1TP in the TP- 
AGB phase. The uncertainties in M c i TP and Mf are propagated to determine o/±m c - The fractional 
contribution of the TP-AGB phase to final core mass (/m c ,tp-agb) is then AM C /Mf. Table [2] lists 
the Mi, Mf, AM C , and /m c .tp-agb of each cluster. 
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2.6. Fuel Consumption: Mass-Light Coupling 



The total en ergy during any evolutio nary phase of a star's life is directly proportional to the 
fuel consumed (IRenzini & Buzzonilll986|) . As nuclear burning is a source of energy for stars, the 
core's growt h during the TP-AGB p hase represents a direct lower bound on the fuel burned during 
this phase (IMarigo & Girardill200l|) . True, there are contributions to the light output from gravita- 
tional contraction of the core and neutrino loses. Ho wever, these can be igno red in our analysis as 
they only represent corrections on the percent level (IMarigo & Girardil 1200 1|) . In addition, signif- 
icant observable light will come from nuclear burning reactions whose products (notably He) are 
expelled in stellar winds and do not end up in the remnant. Hence, we define the growth of the 
core as a strict lower bound on the fuel consumed in the TP-AGB phase. 

Mathematically, w e couple the consumed fuel to the energy released during the TP-AGB 
phase via (cf. eq. 5 from IMarigo & GirardfeOOll) : 



F(M f ,Mi) = -?- / L Mi (t)dt. 



(1) 



The fuel, F(M,), is expressed in solar masses and the integral is the sum of the energy released 
during the TP-AGB phase. The conversion from energy to mass is represented by A H , the efficiency 



of H b urning reactions. Here, we adopt A H = 9.75 x 10 10 L & yr MJ following IMarigo & Girardi 



(120011) . Above, we note that nucleosynthesis is the dominant stellar energy source. The fuel 
consumed in the TP-AGB phase is therefore: 



Ftp-agb ^ (X ia + 0.l)(AM c )+M; F - AOB (He) + a.l-Y')M;(CO). 



(2) 



X12 is the surface mass fraction of hydrogen after the second dredge-up event, Y' is defined as 
^1,2/(^1,2 + ^1,2) where Y\ 2 is the mass fraction of helium in the envelope after the second dredge- 
up, and My P ~ AGB (He) and My(CO) refer to the stellar yield of He and CO, respectively, during 
the TP-AGB phase. We assume an initial H abundance (X q) of 0.71 and map to X lt2 using the 
tabulated surface compositions in lBertelli et al.l (I2008L |2009|) . Irrespective of our lack of knowledge 
regarding the yields of He and CO we can construct our strict lower limit using only the observed 
growth of the core during the TP-AGB phase. The minimum bound on the fuel consumed during 
this phase, and hence the energy released, is obtained by setting the fuel in equation [2] equal to the 
first term on the right hand side and setting that equal to F(M{) in equationQ] Solving for J L Mi (t)dt 
in equation [H we obtain the minimum energy release (£ m i n ) necessary to increase the core mass 
by AM C . The uncertainty in £ min is propagated from a Fj 
Tabled 



, . We list this value for each cluster in 



We compare £ min to th e predicted integrated luminosity of evolutionary tracks incorporat- 
ing the TP-AGB phase (e.g. [pietrinferni et al.l 120041 : IMarigo & Girardil I2OO7I : iBertelli et al]l2008[ 
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20091). Straightforwardly, we can integrate J L Mi (t)dt in the TP-AGB phase for any choice of 
model and find the predicted energy release. As they incorporate up to date physics and produce 
tracks that span a wid e range of metallicity, we use the latest TP-AGB models from the PADOVA 



group in our analysis ([Bertelli et al.1 12008. 2009). The range of composition in our cluster sample 
is covered by their scaled-solar tracks with (Z=0.008, Y=0.26), (Z=0.017, Y=0.26), and (Z=0.40, 
Y=0.30). We create a grid of J L(t)dt in the TP-AGB phase as a function of M Y and composi- 
tion. We interpolate the grid quadratically in mass and lineally in metallicity to determine the total 
predicted luminosity output ( J L M . [Fe / H ^(t)dt). The procedure to calculate the central value for 
/ ^Mi,[Fe/H](0^ an d c/l m . We / m (t)dt is similar to the calculation of Mi and its error (Section [273b . The 
ratio f(J L Mi ,[Fe/H](t)dt) = £mi n / f L Ml ,[Fe/H](t)dt represents the fraction of predicted energy released 
in the TP-AGB phase accounted for by core growth. The quantity 1 — /(/ L M . [Fe / H ](0^0 is now 
the fraction of the total predicted energy that is traced by the products of relic Helium and Carbon 
burning released to the ISM during the TP-AGB phase. Given either theoretical or observational 
constraints on the total light output in the TP-AGB phase, our analysis leads to direct predictions 
of TP-AGB stellar yields. 

In the above prescription we neglect neutrino loses and the the contribution of the core's grav- 
itational contraction to the energy budget of the TP-AGB phase. Gravitational contraction can 
provide, at most, ~ 5% of the energy budget of the TP-AGB phase. Nucleosynthesis of H into 
He and eventually C and O has an efficiency of ~ 5 x 10 18 erg jg whereas the gravitational energy 
released as luminous radiation per unit mass from dwarf creation is ~ 2 x 10 17 erg/ g. Energy gen- 
eration from core contraction is negligible given the uncertainties in the dominant energy source- 
nucleosynthesis. We determine the energy output accounted for by the growth of the core and 
establish constraints on the stellar yields of He, C, and O during the TP-AGB phase as a function 
of initial stellar mass in Sections [331 and 13.41 respectively. 



3. Results 

Observations of white dwarfs in open clusters combined with established evolution models 
constrain the fuel consumption, and hence light output, during the TP-AGB phase. Our results 
demonstrate that the stellar core grows by a non-negligible amount during the TP-AGB phase, 
regardless of progenitor mass. Below, we address the evolution of core mass growth in TP-AGB 
stars as a function of initial mass, couple this core mass increase to a lower bound on TP-AGB star 
light output, compare our results with recent theoretical models of the TP-AGB phase, and discuss 
how our results change with different overshooting conditions. 
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3.1. Initial Final Mass Relation 



We plot the IFMR for our sample in Figure Q] The observed, final masses and their calcu- 
lated progenitor masses are represented by the gray points. The weighted mean Mf-M\ pair for 
each cluster is in black. When compared with the literature from which the initial data came, 
our final masses are within the errors of earlier measurements. Typically, the final masses of a 
given cluster are consistent across different IFMR studies despite coming from various theory 
groups (e.g. IWilliams et al.l l2009a|) . On the other hand, the cluster parameters that dictate the 
calculated initial masses are less well constrained; cluster age, in particular, is notoriously diffi- 
cult to determine. The uncertainties in the inputs to MS isochrone fitting: namely composition, 
reddening, and distance, all contribute to uncertainty in the cluster age and hence, progenitor 
masses. Progenitor masses calculated in our analysis rather than in S09 still show good agree- 
ment with those of other studies assuming the same cluster ages, e.g ., in M37 our weighted mean 
Mj= 3.21 ±0.1 1M for r clus = 563 ± 3 1 Myr and lFerrario et al]J2005n . who use r clus = 520 ± 80 Myr, 
find Mj= 3.27 ±O.12M . Irrespective of the assumption of convective overshooting, the average 
fractional uncertainty in our sample for M f (Mj) is 3% (7%). 

The red line in Figure Q] connects the anticipated core mass at the first thermal pulse (M c i T p) 
fo r each cluster's weighte d mean initial mass. We determine M c j T p given M; and [Fe/H] according 
to IPietrinferni et al.l (|2004f) (see Section [23]). The core mass at the first thermal pulse is a mono- 



tonically increasing function of initial mass for a given metallicity. The clusters in our sample span 
0.32 dex in [Fe/H]. The differing compositions of the two clusters with the lowest progenitor mass 
explains the small dip in the red line at Mj~ 2M Q . The relationship between the core mass at the 
first thermal pulse and the progen itor mass is a robus t theoretical prediction; comparison between 
the models used here and those of iGirardi et al.l (|2000|) yield identical results to within 2%. In Fig- 
ured] all the cluster mean final masses are larger than the predicted core mass at the first thermal 
pulse. Given nominal assumptions regarding stellar evolution, current cluster observations demand 
that the core grows substantially during the TP-AGB phase; thus, TP-AGB stars must contribute 
significantly to cluste r luminosity. While this idea has been gaining traction in the literature (e.g. 
Maraston et al.ll2006|) . we show, using the IFMR and accepted stellar interior theory, that TP-AGB 
stars must be considered a large luminosity source in population synthesis models. 



3.2. Core Growth in the TP-AGB Phase 

The core growth during the TP-AGB phase is a strict lower limit on the mass of nuclear 
burning products in TP-AGB stars and the requisite contribution to stellar populations' bolometric 
luminosity. In Figure [2] we illustrate the fraction of the final core mass grown during the TP- 
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AGB phase as a function of initial mass. We define this fractional mass as /m c ,tp-agb= AM C /M t - 
where AM C is the difference between the final remnant mass and the core mass at the onset of 
the TP-AGB phase. The data points are made using the weighted mean initial and final masses 
for each cluster. As mentioned in Section [2] we compute the progenitor mass and M c 1TP using 
evolutionary tracks both with and without overshooting. Filled squares represent masses found 
with models and cluster ages incorporating overshoot; open squares are masses that do not take 
overshooting into account. Dashed lines connect the results of these two cases for each cluster. The 
error bars, calculated according to Sections 12.51 and 12.31 represent one sigma errors in fractional 
core mass gained during the TP-AGB phase and initial mass, respectively. The red lines connect 
the moving weighted mean of each data point type, illustrating the overall trends in the figure. 
At every filled or open point, we calculate the weighted mean of itself and its closest neighbor in 
both directions of Mi (end points are only averaged with the closest data point). The solid red line 
connects the averages of the filled data points (OS); the dashed red line connects the open square 
(nOS) averages. 

These moving averages highlight the large increase in fractional core mass growth during 
the TP-AGB phase between 2 and ~ 3M Q (up to 4M in the no overshooting case). In the case 
where convective overshooting is considered (OS), the weighted mean /m c ,tp-agb= 0.20 ±0.01 
when 1.9 < Mj < 3.6 and drops by a factor of two to 0.08 ±0.02 elsewhere. This broad peak 
in /m c ,tp-agb above the remainder of the sample is significant at the 3.5a level. For progenitor 
masses calculated with models and cluster ages that do not include convective overshooting (nOS), 
Aic, tp-agb is generally smaller than its OS counterpart and the fraction of core mass gained in the 
TP-AGB phase peaks at a slightly higher progenitor mass (3.5M Q versus 3M Q ). These shifts are a 
natural consequence of the inclusion or absence of convective overshooting in evolutionary models 
and isochrones. Convective overshooting models predict a larger He core and longer lifetime 
along the main sequence compared to nOS models. Due to these larger cores, M c 1TP is larger 
for a given M x when overshooting is considered. If the initial masses calculated with and without 
convective overshooting were the same, /m c ,tp-agb would be larger in the nOS case. However, OS 
isochrones assign older cluster ages (and subsequently smaller Mi) than nOS isochrones given the 
same cluster CMD. The M c iT p - Mi relation is monotonically increasing and thus the nOS M c j T p is 
higher than its OS counterpart in the same cluster when Mj< 5M Q (the M c 1TP - Mi relation flattens 
slightly at higher M ; ; nOS M c 1TP are similar to the OS M c 1TP in this mass range). Still, the nOS 
average curve shows the same shape as that composed of OS progenitor masses. On a cluster by 
cluster basis, the OS and nOS core growth fractions are typically within 5% of each other. While 
observations of gaps in open cluster CMDs and the HR diagram postions of binary stars give 
more credance to models invoking core overshoot, the precise extent of convective cores is still a 
matter of active debate. It is thus encouraging that our results are surprisingly independent of this 
significant theoretical uncertainty. Given the similarity of our results regardless of overshooting 
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behavior and the facts that models with convective overshooting are favored, future sections will 
only address our OS results. 

Interestingly, the rise in /m c ,tp-agb r oughly coincides with th e peak in the lifetime function 



predicted by recent TP-AGB models (e.g. Marigo & Girardil 120071) . The lifetime of the TP-AGB 



phase is set by the time between the first thermal pulse and the removal of the envelope. The 
envelope can be removed via mass loss to the ISM or through outward advancement of the H- 
burning shell (hot bottom burning). In many TP-AGB models, the lifetime of the phase depends 
critically on the occurrence of the third dredge-up (3DU) episode. In the 3DU, carbon is mixed into 
the envelope, causing the C/O ratio to rise and the opacity of the envelope to increase increase. 
In turn this creat es and sustains a dust-dr iven wind that greatly enhances the mass loss rate of 



the envelope. In iMarigo & Girardil (|2007[) . the 3DU occurs when the mass of the core exceeds 
a threshold (M crit ) that, at solar metallicity, is relatively flat when M ; < 3M and becomes a steep 
function of M\ at higher masses. Though A/ c j T p typically grows with Mi, the M c i TP -M; relationship 
plateaus between 1.5M Q and 3M Q . The core must grow fractionally more in the TP-AGB phase 
to reach M crit in this mass regime; therefore, TP-AGB star lifetimes are longest for these initial 
masses. More massive stars (i.e. > 3M ) experience their third dredge-up more quickly as their 
core mass at the onset of the TP-AGB phase is closer to M crit . Conversely, we observe that when 
Mj<~ 1.8M Q , stellar cores at the first thermal pulse are already a substantial fraction of M crit . The 
lowest Mi in our sample, 1.75M Q , has the smallest /m c ,tp-agb and is predicted to have a shorter 
lifetime than 2 M Q and 3 M TP-AGB stars. 



3.3. Fuel from Stellar Core Growth 

Using the procedure outlined in Section 12.61 we determine the fuel consumption during the 
TP-AGB phase required by core mass growth. We have already determined the mass added to 
the core during the TP-AGB phase (AM f ) for each cluster. Using equations [2] and Q3 AM C is 
converted to its equivalent energy output (£ ra i n ). We have shown that the theoretical uncertainty, 
dominated by the chosen efficiency of convective overshoot, associated with the core mass at the 
first thermal pulse is relatively small. Thus, C min is a strong lower bound on the total energy output 
of TP-AGB stars stemming almost entirely from observational constraints. We compute C min as 
a function of initial mass for each cluster (Table [2]). Core growth in TP-AGB stars demonstrates 
that the TP-AGB phase contributes a significant portion of a intermediate mass star's overall light 
generation. Our results emphasize the strong link between the observables of the IFMR and the 
lower limit on the integrated light during the lifetime of TP-AGB stars. 

We now quantify the level of agreement between £ m j n and the predicted total light output 
by current evolutionary models of TP-AGB star. We define /e,am c = -Cmin/ / L Mu [ Fe / H ^(t)dt where 
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/ ^Mi,[Fe/H](0^ is the predicted total energy in the TP-AGB phase and is obtained by integrat- 
ing TP-AGB stellar tracks . We choose the latest PADOVA group TP-AGB tracks, produced by 



Bertelli et al.1 (120081. 120090 . as they make predictions for a wide range of metallicities and incor- 



porate up to date physics. In Tabled we list /e,am c and its error on a cluster by cluster basis. 
The weighted mean of /e.am c in our sample is 0.56 ±0.04. The individual /e,am c for each cluster 
is within 2a of this mean while all but the cluster with the sma l lest progenit or mass, NGC 6819, 
and Praesepe are within la of this mean. Unless iBertelli et al.l (120081 12009) substantially under- 
represent the energy output of the TP-AGB phase, the IFMR provides an estimate of the total light 
emitted in the TP-AGB phase accurate to within a factor of ~ 2. 



3.4. Fuel from Stellar Yields 



If we make the reasonable assumption (Section 12.61) that the available fuel in the star during 
the TP-AGB phase comes primarily from nuclear burning, the quantity 1 — /e,am c is the fraction 
of total predicted energy accounted for by nucl eosynthetic p r oducts that w ere expelled into the 



ISM rather than added to the core. Assuming IBertelli et al.1 (120081. 12009Q correctly predict the 



total energy output in TP-AGB stars, the product (1 -/e,am c ) * / L Ml ,[Fe/H](t)dt is the amount of 
light unaccounted for by core mass growth. Setting (1 -/e.amJ * / L Ml ,[Fe/H](t)dt = / L Mj (t)dt in 
equation [H we obtain the available fuel, in solar masses (i r TP-AGB,He) 5 that comes from the stellar 
yield of the star. In equation [21 we make a distinction between the yield of He and that of CO. In 
practice, the binding energy of He is 1/10 that of CO; thus, me make the approximation that He 
is the sole component of the stellar yield. The calculated Fxp-AGB.He from each cluster is shown 
in Table [2l We find our sample's weighted mean i 7 T p-AGB,He= 0.08 ±0.01 M Q . Every individual 
measurement of the stellar yield is within two a of this mean. These results indicate that every star 
born at ~ 1.5M < M { < 6M Q deposits ~ 0. 1 M & of He into the ISM during the TP-AGB phase. 
We discuss potential ramifications and tests of this prognosis in the Section HI 



4. Summary and Discussion 

Using the WD IFMR and robust stellar evolutionary theory, we have shown that a significant 
mass fraction of the final, stellar core is generated during the TP-AGB phase. Through simple 
fuel consumption arguments, the TP-AGB phase must therefore emit a substantial amount of light- 
proportional to the aforementioned core growth. There is measurable core growth during the TP- 
AGB phase in all nine clusters in our sample. The representative progenitor masses range from 
~ 1.7 < M Q < 6 and they span ~ 0.30 dex in [Fe/H]. The least constrained theoretical assumption 
in stellar evolutionary models up to the onset of the TP-AGB phase is the depth of convective 
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overshooting and whether or not it takes place. Our main result is valid over a broad swath of 
theoretical input parameter space, including the uncertain nature of core overshooting. Therefore, 
we have constructed a strict lower bound on the light output during the TP-AGB phase as a function 
of initial mass that is predominantly dependent on observational constraints from white dwarfs and 
their host cluster properties. 

The fractional contribution of TP-AGB stars to the total light emitted over their lifetime peaks 
when Mi ~ 2-3.5M . In general, ~ 20% of the core is built during the TP-AGB phase in this 
mass range, double that of lower and higher mass stars. The rise and fall of the TP-AGB's impor- 
tance with initial ma ss generally agrees with the relative lifetimes of TP-AGB stars predicted b y 
current models (e.g. birardi & Marigoll2007l : llVIarigo & Girardill2007l : bertelli et allbood. 12009k 
The maximum of the TP-AGB fractional energy output at ~ 2-3.5M Q has strong implications for 
interpreting the integrated light of elliptical galaxies. The Infrared light from TP-AGB stars will be 
a one of stronger principle components of the spectral energy distributions of galaxies harboring 
stellar populations 2 to 3 Gyr old. 

Constraining the importance of the TP-AGB phase via the IFMR becomes increasingly dif- 
ficult for progenitor masses above ~ 5M . High mass stars have short lifetimes, implying that 
Tcius -Tcoo\ is a small number. In that case, nominal errors in r c i us and r coo i have a fractionally larger 
impact on the uncertainty in progenitor age at high mass than at low mass. Additionally, different 
convective overshooting efficiencies produce more discrepant core masses at various stages of evo- 
lution for higher initial masses. Still, we find compelling evidence that ~ 10% of the final core is 
built up during the TP-AGB stage at high mass regardless of convective overshoot efficiency. The 
upper limit in progenit or mass for stars expe riencing the TP-AGB phase is still controversial and 
empirical data is scare (iWilliams et al.ll2009b|) . While this work suggests 5.5 M Q stars experience 
a TP-AGB phase, the IFMR will provide a strong means to determine the highest initial masses to 
become TP-AGB stars if progenitor masses in this regime can be better constrained in the future. 

Our results have potentially important consequences for chemical evolution models. We de- 
termine the total energy output in the TP-AGB phase required by the nuclear reactions that directly 
lead to core mass growth. Assuming nucleosynthesis is the only source of energy generation in 
the TP-AGB phase, the only other tracer of energy output would be nucleosynthetic byproducts 
expelled to the ISM. Any measurement of the total light output of a TP-AGB stellar population, 
combined with this work, places direct constraints on the yield of that population (principally he- 
lium). If the helium yield is significant, it would impact the AF/AZ relationship. Additional 
helium in the ISM from TP-AGB stars wou ld alter the luminance and lifetimes of subsequent 
generations of stars. Assuming the models of iBertelli et al.l (|2008L 12009) correctly predict the to- 
tal light in the TP-AGB phase, we find that 0.08 ±0.01 M Q of helium are released into the ISM 
by TP-AGB stars. TP-AGB stars are short lived; as such, observations of TP-AGB populations 
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and limits on their luminosity have proven difficult in galaxies more distant than the Magellanic 
Clouds. However, the above argument can be reversed: given limiting case assumptions of the 
helium yield of TP-AGB stars and IFMR data, one can establish both upper and lower bounds on 
the integrated light from TP-AGB stars. 

Early evidence suggests that modern models of TP-AGB stars make realistic predictions for 
their total light output. Mass addition to the core accounts for ~ 5 0% of the integrated luminosity 
from the TP-AGB phase predicted by iBertelli et all J20081 l2009h . If / E AMc had exceeded 100%, 



we would conclude that the chosen models underestimated the amount of light emitted by TP- 
AGB stars. Alternatively, if /E,AM c was small (< 10%), we would calculate a relatively large stellar 
yield. In the future, empirical evidence describing the behavior of the AY / AZ relationship or 
integrated light measurements of TP-AGB populations will constrain TP-AGB stell ar yields and 



impose tight er restrictions on the upper bound of total energy release. Currently, the Bertelli et al 



(120081 12009) models predict a n amount of luminous ene rgy release by TP-AGB stars that is in 



agreement with observations. Bertelli et al. ( 2008 . 2009|) predict that, at its peak, the TP-AGB 



phase is responsible for ~ 33% of the total light output of stars with 2.2 M e < Mj < 2.5M ; this 
fractional contribution falls to 20% when M { = 3.5 M Q and is ~ 10% at M { = 4.0 M Q . Though our 
current results and errors cannot rule it out, it is unlikely that TP-AGB stars are responsible for as 
much as 40% of cluster light unless the stellar yields of TP-AGB stars are significantly higher than 
those predicted in Section [3~4l 

This work suggests additional tests of population synthesis models. The predicted light emis- 
sion from stellar populations in SPS models provides a upper bound on the expected remnant mass 
as a function of initial mass. In the future, population synthesis codes should quantify their agree- 
ment with the IFMR. Various theoretical question marks in SPS modeling can be answered, or 
at least constrained, by the IFMR. The degree of convective overshooting in stellar evolutionary 
models has an appreciable impact on implied progenitor masses and the M c 1TP -Mi relationship; 
hence, for a predicted light output in the TP-AGB phase, there will be a statistically significant 
difference between the remnant mass functions predicted with different overshooting parameters. 
If we demand that M c j TP is such that the core grows significantly during the TP-AGB phase, we 
can place instructive bounds on convective overshoot prescriptions for some clusters. Additionally, 



the l uminosity of TP-AGB populations represent a major uncertainty in SPS models ( Conroy et al 



2009). Our results showing core mass growth in the TP-AGB phase as a function of initial mass 
represent a floor to TP-AGB star luminosity and can already exclude some limiting case models. 
Obviously, models assuming that TP-AGB stars do not contribute to a population's integrated light 
can no longer be physically motivated. As mentioned above, this work and assumptions as to the 
chemical yields of TP-AGB stars is sufficient to constrain TP-AGB light output. These new con- 
straints are valuable to future SPS models and will greatly reduce the uncertainties associated with 
the one of the most theoretically dubious stellar evolutionary phases.. 
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The most prominent sources of error in this analysis are the WD and cluster observations. 
Measurements of WD surface gravities and temperatures provide constraints on remnant masses 
and cooling times while cluster observations yield the cluster age and composition. Uncertainties 
in theory linking these observations and the parameters of interest are typically small compared 
to the observational errors (S09). By assuming either perfect WD or cluster measurements and 
repeating our analysis, we find that uncertainties in these two types of observables contribute sim- 
ilarly to the resulting error in core mass growth during the TP-AGB phase. More precise WD or 
cluster measurements will provide stricter bounds on the energy output of TP-AGB stars. Another 
potential method to improve the precision of the IFMR (and our results) would be a survey specif- 
ically designed to find young white dwarfs in open clusters. In any given cluster, the youngest 
white dwarfs would correspond to the oldest progenitor lifetimes, reducing the impact of uncer- 
tainties in the WD and cluster observations on the implied initial mass and M c 1TP . More precise 
measurements of the IFMR would directly lead to more precise calculations in our procedure. 

Determining the metallicity dependence of TP-AGB star energy output is an intriguing next 
step. However, there are not enough open clusters at extreme metallicities to empirically calibrate 
such a dependence. We would need to increase the number of clusters in our sample, even if the 
overall spread in metallicity does not change. With many more clusters, one could reproduce our 
results in bins of metallicity. After empirically characterizing how the energy output of the TP- 
AGB phase scales with composition (albeit over a limited range), theory may be able to provide 
a physical model for TP-AGB star luminosity as a function of initial mass and stellar density. If 
said model could be extrapolated to low metallicity, it would provide a theoretical constraint on the 
contribution of TP-AGB stars to the spectral energy distributions of high redshift galaxies. At high 
metallicity, the lowest mass stars may go dir ectly to the white dwarf stage, missing the core He 



burning stages and beyond iKilic et al.ll2007l) . This effect is not currently included in population 
synthesis models and might be important for the low redshift properties of giant elliptical galaxies. 
Physically motivated priors on TP-AGB star light output would be a substantial advancement of 
semi-analytic galaxy evolution models as SPS codes still have widely varied characterizations of 
these stars. 

We have shown that the WD IFMR places powerful constraints on the energy released during 
the TP-AGB stage of stellar evolution. Future observations of TP-AGB populations will constrain 
the helium yields of these stars- an important factor for chemical evolution models to consider. Al- 
ternatively, nominal assumptions regarding the yield of these stars, in conjunction with this work, 
result in a narrower range of possible TP-AGB population luminosity than previously considered. 



We would like to thank David Weinberg for his encouragement and thought-provoking dis- 
cussion. 
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Table 1 . Cluster Metallicity, Distance, and Age 



Cluster 


[Fe/H] 


(m-M)o (mag) 


OS Age (Myr) 


nOS Age (Myr) 


References 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 6819 


0.09 ±0.03 


12.15±0.20 


2000 ±200 


1500 ±200 


1,2 


NGC 7789 


-0.04 ±0.05 


12.19±0.12 


1442 ±100 


1042 ± 100 


3,2 


Hyades 


0.13 ±0.06 


3.33±0.01 


640 ±40 


440 ±40 


4,2 


Praesepe 


0.14 ±0.02 


6.33 ±0.04 


637 ±50 


440 ±40 


5 


M37 


0.04 ±0.04 


10.83±0.13 


554 ±54 


354 ±43 


6 


NGC 3532 


0.02 ±0.06 


8.48 ±0.20 


316±80 


216±80 


4,7 


NGC 1039 


0.07 ±0.04 


8.58±0.15 


250 ±25 


150 ±30 


8,2 


NGC 2516 


-0.04 ±0.05 


8.03 ±0.04 


137 ±29 


91 ±26 


5 


M35 


-0.21 ±0.10 


10.48 ±0.12 


124 ±30 


88 ±25 


9,2 



Note. — Cluster compositions, distances, and ages. (1): Cluster. (2): [Fe/H]. (3): Absolute 
distance modulus. (4): Cluster age using isochrones with overshooting. (5): Cluster age using 
isochrones without overshooting. (6): Reference for cluster metallicity and distance; if two refer- 
ences are listed, the first applies to metallicity and the second, distance. Note that the distances in 
column 3 may be modified from their references if different cluster compositions were adopted. 



References. 



Ref erenc es: 1 . iBragaglia et al 



Tautvaisiene et al. (2005); 4. 



Gratton 



( 


2001 


); 2. 


S alar is et al. 



2009b; 

(120001): 5. lAn et al.l d2007l): 6. lHartman et al.l (L2008); 



3. 

7. iKharchenko et all d2005h : 8. ISchuler et all d2003h : 9. iBarrado y Navascues et all d200ll) 
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Fig. 1. — The IFMR for the 48 WDs in our sample. The progenitor masses are calculated using 
cluster ages and stellar evolutionary theory that incorporate a degree of convective overshooting 
(IPietrinferni et al.ll2004|) . The M r M x pairs for each WD are in light gray; error bars represent one 
a deviations in each quantity. The open black squares show the weighted mean Mf and M\ for the 
nine clusters in our sample. Error bars are the one a error of the means. Given an initial mass and 
cluster metallicity, theory predicts the mass of the core at the first thermal pulse (M c j TP ). The solid 
red line connects the nine M c 1TP - M[ points created using the weighted mean Mj and the [Fe/H] 
of each cluster. The observed remnant masses in these clusters all lie above the red line. As M f is 
the mass of the core at the tip of the AGB, stellar cores must grow during the TP-AGB phase. 
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Fig. 2. — The fraction of the final remnant mass built up in the TP-AGB phase as a function of 
initial mass. The squares are calculated using the weighted mean initial and final masses of each 
cluster. Filled squares use evolution theory and isochrones incorporating convective overshooting; 
open squares represent the limiting case of no overshooting. The dashed lines connect these two 
choices for each cluster. The solid (dashed) red line connects the results of a moving weighted 
mean of the filled (open) squares. The fractional core mass growth in the TP-AGB phase has a 
broad peak between 2-3M . 



-25- 



Table 2. Fractional Core Growth, Corresponding Energy Output, and Estimated Helium 

Yields 



Cluster 


Mj (M ) 


M f (M ) 


fc,TP-AGB 


AM C (M ) 


£min (1O 9 L JT) 


/e.AM c 


^TP-AGB.He ( M Q ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC6819 


1.75 ±0.08 


0.55 ±0.02 


0.06 ±0.03 


0.03 ± 0.02 


2.66 ± 1.22 


0.29 ±0.14 


0.07 ±0.01 


NGC7789 


1.94 ±0.05 


0.61 ±0.02 


0.18 ±0.04 


0.11 ±0.02 


8.29 ± 1.65 


0.56 ±0.14 


0.07 ± 0.02 


Hyades 


2.86 ± 0.03 


0.71 ±0.01 


0.20 ±0.02 


0.15 ±0.01 


1 1.07 ± 1.01 


0.56 ±0.05 


0.09 ±0.01 


M37 


3.10 ±0.07 


0.80 ±0.03 


0.23 ± 0.04 


0.18 ±0.03 


I3.7S ± 2.21 


0.71 ±0.12 


0.06 ± 0.02 


Praesepe 


3.21 ±0.16 


0.79 ± 0.02 


0.20 ±0.05 


0.15 ±0.04 


11.70 ±2.74 


0.62 ±0.14 


0.07 ± 0.03 


NGC 3532 


3.57 ±0.21 


0.85 ±0.02 


0.16±0.05 


0.13 ±0.05 


10.07 ±3.40 


0.56 ±0.19 


0.08 ± 0.04 


NGC 1039 


4.00 ± 0.23 


0.87 ± 0.04 


0.09 ± 0.06 


0.08 ± 0.05 


5.97 ±3.85 


0.41 ±0.27 


0.09 ± 0.04 


M35 


5-14 ± 0.31 


0.96 ±0.03 


0.10 ±0.04 


0.09 ± 0.03 


6.42 ± 2.42 


0.70 ±0.43 


0.03 ± 0.04 


NGC 2516 


5.47 ±0.59 


0.99 ± 0.03 


0.13 ±0.06 


0.12 ±0.06 


8.59 ±4.37 


0.96 ±0.66 


0.00 ± 0.06 


Weighted mean 


2.65 ± 0.02 


0.73 ±0.01 


0.16±0.01 


0. 12 ± 0.01 


8.37 ±0.60 


0.56 ± 0.04 


0.08 ±0.01 



Note. — Initial mass, final mass, fractional core growth and corresponding energy output, and estimated helium yields of each cluster. (1): Cluster. (2): Progenitor Mass in 
Mq . (3): Remnant mass in Mq . (4): Fraction of remnant mass added to the core in the TP-AGB phase. (5): Core mass added in the TP-AGB phase in Mq . (6): Minimum total 
energy output required by TP-AGB core growth in 1 Q 9 L^yr. (7): Fraction of e nergy in TP-AGB phase predicted by Bertelli et al. i 2008, 2009) accounted for by core growth. (8): 
Assuming the total energy output of TP-AGB stars in lBertelli et ai]|2008ll2009l . estimated stellar yields of helium in Mq . 



